Theory of the striped superconductor 
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We define a distinct pliase of matter, a pair density wave (PDW), in wliich the superconducting 
order parameter, (p{r,f'), varies periodically as a function of position such that when averaged over 
the center of mass position, (r + f")/2, all components of (j) vanish identically. Specifically, we study 
the simplest, unidirectional PDW, the "striped superconductor," which we argue may be at the 
heart of a number of spectacular experimental anomalies that have been observed in the failed high 
temperature superconductor, La2-a:Baa;Cu04. We present a solvable microscopic model with strong 
electron-electron interactions which supports a PDW groundstate. We also discuss, at the level of 
Landau theory, the nature of the coupling between the PDW and other order parameters, and the 
origins and some consequences of the unusual sensitivity of this state to quenched disorder. 



I. A NEW PHASE OF MATTER 

Superconductivity which arises from the pairing of 
electrons in time reversed states is well understood 
as a weak coupling "Fermi surface" instability - a 
consequence of arbitrarily weak effective attractive 
interactions fi In this paper we explore a new type of 
superconducting state, a "pair density wave" (PDW), 
which is a distinct state of matter, and which does not 
occur under generic circumstances in the weak coupling 
limit - it requires interactions in excess of a critical 
strength. The PDW spontaneously breaks the global 
gauge symmetry, in precisely the same way as a conven- 
tional superconductor does. Thus, the order parameter 
is a charge 2e complex scalar field, (j). However, it also 
spontaneously breaks some of the translational and point 
group symmetries of the host crystal, in the same way as 
a conventional charge-density wave (CDW). In a PDW, 
for fixed f—f', the order parameter 



A(r,r') = (^|(r)^|(r")) 



(1.1) 



is a periodic function of the center of mass coordinate, 
R = {f+r')/2. Here, 'ipl-{r) is the electron creation 
operator at position r with spin polarization a. We will 
focus on the case in which translational symmetry is bro- 
ken only in one direction, i.e. a unidirectional PDW or 
equivalently a "striped superconductor." 

Accompanying the PDW there is, as we will show be- 
low, induced CDW order with half the period of the 
PDW. However, the PDW differs from a state of coex- 
isting superconducting and CDW order—, which has also 
been previously called a "pair-density wave" ;^ Unlike the 
pair-density- wave state of Ref. [3] , the average value of the 
superconducting order parameter vanishes for the PDW 
state we discuss here. This is a defining symmetry prop- 
erty of this state. The PDW is more closely analogous 
to the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state^ 
which arises, under appropriate circumstances, when the 
electron gas is partially polarized by an applied magnetic 
field. However, explicit time-reversal symmetry break- 
ing is an essential ingredient of the FFLO state, and 



is responsible for lifting the degeneracy ("nesting") be- 
tween time-reversed pairs of quasiparticle states. In the 
absence of quenched disorder, the PDW state preserves 
time-reversal symmetry.^ 

In general, there is no necessary relation between spin- 
density wave (SDW) and PDW order - whatever rela- 
tion there is derives from common microscopic physics 
rather than from general symmetry conditions. How- 
ever, the two orders can be linked if one postulates a 
larger, emergent S0(5) symmetryi^ which unifies the 
superconducting and SDW orders. In this context, an 
early speculation^^ concerning the existence of a unidi- 
rectional SO (5) spiral density wave state is particularly 
interesting.^"* Such a state would consist of interleaving 
SDW and superconducting spirals, both with the same 
period. The fact that this state is a spiral implies that 
the phase of (j) varies as a function of position, and hence 
that there are equilibrium currents. Thus, the S0(5) 
spiral, though similar to a state of coexisting PDW and 
coUinear SDW order, is not the same. We have not found 
any microscopic model that supports a spiral-PDW phase 
- whether coexisting with SDW order, or not. However, 
it is an interesting state and deserves further studyJ^ 
Many features of the superconducting state depend only 
on the existence of a charge 2e order parameter. This 
includes zero resistance, the Meissner effect, flux quanti- 
zation and, in 3D, the existence of a finite critical current. 
Therefore, these properties are expected in a PDW (in 
the absence of quenched disorder). 

However, other properties are very different. For in- 
stance, in an s-wave superconductor, there is a full gap 
in the quasiparticle spectrum. For a d-wave superconduc- 
tor, there are typically nodal points on the Fermi surface 
which support gapless quasiparticle excitations, but the 
density of states still vanishes at the Fermi energy, even 
in the presence of coexisting CDW or most forms of SDW 
order*^. In contrast, in a PDW portions of the Fermi sur- 
face are typically ungapped, which as in the case of many 
CDWs, results in the reconstruction, but incomplete gap- 
ping of the underlying Fermi surface^^J*^. In a conven- 
tional superconductor, the anisotropy of the supercon- 
ducting coherence length is determined by the anisotropy 
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of the single-electron effective mass; in a PDW, depend- 
ing on the precise character of the translation symmetry 
breaking, the anisotropy of the superconducting state can 
be parametrically larger than the normal state resistivity 
anisotropy. 

Perhaps the most significant new feature of a PDW 
superconductor is its anomalous sensitivity to (non- 
magnetic) quenched disorder, in stark comparison to a 
conventional superconductor. Specifically, although the 
disorder potential does not couple directly to the super- 
conducting order parameter, it produces random varia- 
tions in the magnitude and sign of the local superfluid 
density. Hence quenched disorder can readily drive a 
PDW into a superconducting XY glass phase. (A similar 
effect was recently predicted^ in a d-wave superconduc- 
tor near the critical point of the superconductor to metal 
transition.) 

Table [J characterizes the PDW, CDW, and uniform 
superconducting (SC) phases by their order parameters 
and their sensitivity to quenched disorder. Here, Aq and 
Aq are, respectively, the uniform SC and PDW (mod- 
ulated SC) order parameters (where Q is the ordering 
wavevector), and pg and p2Q are the two relevant CDW 
order parameters. By definition, the pure PDW phase 
has non-zero Aq and a vanishing Aq. It is thus different 
from a phase of coexisting SC and CDW order, which has 
a non-zero Aq and Aq (as well as a subsidiary pg order). 
(In various places in Refs. 0, as an aid to intuition, the 
SC-I-CDW and/or the CDW states have also been re- 
ferred to as "Cooper pair crystals") As mentioned above, 
the main new characteristic of the PDW (as opposed to a 
uniform SC) is its sensitivity to disorder. Even weak dis- 
order destroys the long range PDW order, turning it into 
a glassy state. The properties of these phases and their 
inter-relations will be derived from an order parameter 
Landau theory in Section [V] 

In a recent publications^ we argued that the remark- 
able properties of La2-2;Baa;Cu04 near a; = 1/8 filling 
can be explained as a consequence of the symmetries of a 
striped superconducting state. The main purpose of the 
present paper is to develop the theory of this state. 

This paper is organized as follows. In Section [TT] we 
give an overview of the physics the striped superconduc- 
tor La2-a;Baa;Cu04 evidenced in recent experiments^i. 
In Section IIIII we discuss microscopic mechanisms that 
give rise to a PDW. In Section IIVI we present several 
solvable models that shed light on the physical origin of 
a PDW superconducting state. Here we investigate the 
role of different model junctions separating two super- 
conducting regions, including an empty barrier (a trivial 
insulator) and an antiferromagnetic insulating barrier. 
In both cases, under certain conditions, we find that a 
TT phase shift between phases of the two superconduc- 
tors is favored. We then build on these insights to con- 
struct a model with a striped superconducting ground 
state. We discuss qualitatively the quasiparticle spec- 
trum of a PDW state. In Section |V] we present a Landau 
theory of this state and discuss in detail the symmetry- 



TABLE L Summary of the different phases discussed in the 
text and their order parameters, as well as their sensitivity 
to quenched non-magnetic disorder. A / specifies that in a 
particular phase, the corresponding order parameter is non- 
zero. See text for details. 
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dictated couplings between the PDW order parameter 
and the nematic, charge density wave (CDW) and spin 
density wave (SDW) order parameters characteristic of 
a stripe state. In a recent publication, Agterberg and 
Tsunetsugu* discussed a Landau theory for a pair den- 
sity wave which complements the results presented in this 
Section. In Section lVTl we present a statistical mechanical 
model that describes a layer decoupled striped supercon- 
ductor, which is proposed to describe the phenomenology 
of La2-KBaa;Cu04. In Scction lVIII wc discuss the implica- 
tions of this theory. Further details of the Landau theory 
are presented in Appendix [X] The case of coexisting uni- 
form and striped d-wave superconductivity is discussed 
Appendix |B] 



II. STRIPED SUPERCONDUCTOR IN 

LA2-^BA^CU04 

La2-2:Baa;Cu04 is currently the most promising can- 
didate experimental system as a realization of a striped 
superconductor—!^. Considerable indirect evidence in 
favor of the existence of PDW order in this material has 
recently been gathered and compiled in Ref. [23]. The 
putative PDW order has the same period as the unidi- 
rectional (spin-stripe) SDW order which is known, from 
neutron scattering studies, to exist in this materia P"^'^^ . 
Strong evidence that a striped SC can be induced in un- 
derdoped La2_2;Sr2;Cu04 by the application of a trans- 
verse magnetic field, which is also known to induce spin 
stripe order, has also recently been presented in Ref. [26| . 

The behavior of La2_a;Ba2;Cu04 is very striking, and 
rather complex. We will not elaborate upon it here (see 
Ref. [1^). However, there are two qualitative features 
of the data on which we would like to focus: 1) With 
the onset of stripe spin order at 42 K,31 there is a large 
(in magnitude) and strongly temperature dependent en- 
hancement of the anisotropy of the resistivity and other 
properties, such that below 42 K the in- plane charge 
dynamics resembles those of a superconductor, while in 
the c-direction the system remains poorly metallic. The 
most extreme illustration of this occurs in the temper- 
ature range 10 K < T < 16 K, in which the in-plane 
direction resistivity is immeasurably small, while the c- 
axis resistivity is in the 1-10 mil range, so the resistivity 
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FIG. 1: Schematic model of a striped superconductor consist- 
ing of an array of superconducting regions (S.C.) separated 
by correlated insulator regions (I). The Josephson coupling 
between neighboring superconducting regions is J. If J < 0, 
the system forms a striped superconducting (or PDW) state. 

anisotropy ratio is consistent with infinity. 2) Despite 
the fact that many clear signatures of superconductivity 
onset at temperatures in excess of 40 K, and that angle 
resolved photoemission has inferred a "gap" ^^'^° of order 
20 meV, the fully superconducting state {i.e. the Meiss- 
ner effect and zero resistance in all directions) only occurs 
below a critical temperature of 4 K. It is very difficult to 
imagine a scenario in which a strong conventional su- 
perconducting order develops locally on such high scales, 
but fully orders only at such low temperatures in a sys- 
tem that is three dimensional, non-granular in structure, 
and not subjected to an external magnetic field. 

We will see that both these qualitative features of the 
problem are natural consequences of the assumed exis- 
tence of a PDW. Indeed, analogous features have long 
been known to be a feature of the spin ordering in the 
same family of materials^^. Specifically, unidirectional 
spin-stripe order is observed to occur under a number 
of circumstances in the 214 family of cuprate supercon- 
ductors. However, when it occurs: 1) Although the in- 
plane correlation length can be very long, in the range 
100 — 400A, the inter plane correlation length is never 
more than a few A, a degree of anisotropy that can- 
not be reasonably explained simply'^^ on the basis of the 
anisotropy in the magnitude of the exchange couplings. 
2) Despite the presence of long correlation lengths, true 
long-range spin-stripe order has never been reported. It 
will be made clear, below, that interlayer decoupling due 
to the geometry of the stripe order, and the rounding of 
the ordering transitions by quenched disorder can be un- 
derstood as arising from closely analogous considerations 
applied to incommensurate SDW and PDW order—. 



III. MICROSCOPIC CONSIDERATIONS 

At first, the notion that a PDW phase could be sta- 
ble sounds absurd. Intuitively, the superconducting state 
can be thought of as the condensed state of charge 2e 
bosons. However, in the absence of magnetic fields, the 
ground-state of a bosonic fluid is always node-less, inde- 
pendent of the strength of the interactions, and therefore 
cannot support a state in which the superconducting or- 
der parameter changes sign. Thus, for a PDW state 



to arise, microscopic physics at scales less than or of or- 
der the pair-size, ^o, must be essential. This physics re- 
flects an essential difference between superfluids of paired 
fermions and preformed bosons^. 

Our goal in this section is to shed some light on the 
mechanism by which strongly interacting electrons can 
form a superconducting ground-state with alternating 
signs of the order parameter. We will consider the case of 
a unidirectional (striped) superconductor, but the same 
considerations apply to more general forms of PDW or- 
der. We will not discuss the origin of the pairing which 
leads to superconductivity. Likewise, we will not focus 
on the mechanism of translation symmetry breaking by 
the density wave, as that is similar to the physics of 
CDW and SDW formation. Our focus is on the sign 
alternation of (/>. Thus, in much of this discussion, we 
will adopt the model shown schematically in Fig. [1] in 
which we have alternating stripes of superconductor and 
correlated insulator. The system looks like an array of 
extended superconductor-insulator-superconductor (SIS) 
junctions, and we will primarily be concerned with com- 
puting the Josephson coupling across the insulating bar- 
riers. If the effective Josephson coupling is positive, then 
a uniform phase (normal) superconducting state is fa- 
vored, but if the coupling is negative (favoring a tt junc- 
tion), then a striped superconducting phase is found. 



A. Previous results 

So long as time reversal symmetry is neither sponta- 
neously nor explicitly broken, the Josephson coupling, 
J between two superconductors must be real. If it is 
positive, as is the usual case (for reasons that will be 
expanded upon in later subsections), the energy is min- 
imized by the state in which the phase difference across 
the junction is 0; if it is negative, a phase difference of tt 
is preferred, leading to a "tt junction." tt junctions have 
been shown, both theoretically and experimentally, to oc- 
cur for two distinct reasons: they can be a consequence of 
strong correlation effects^i^^i^ or magnetic orderin g^^i'^^ 
in the junction region between two superconductors, or 
due to the internal structure ( e.g. d-wave symmetry) 
of the superconductors, themselves^^-. 

In the present paper, we will build on the flrst set of 
ideas. Until now, tt junctions have been confined to sys- 
tems in which the Josephson tunneling is dominated by 
a single impurity site or quantum dot. It is non-trivial 
to extend this mechanism to the situation in which there 
is an extended barrier in which J is proportional to the 
cross-sectional "area" of the junction^^. Indeed, as we 
shall see below, under most circumstances, even in the 
case in which the Josephson tunneling through an iso- 
lated impurity would produce a negative J, for a barrier 
consisting of an extended area of identical impurities, J 
will typically be positive. One thing that we achieve be- 
low is to obtain a proof that there are circumstances in 
which J of such an extended barrier is negative, and to 
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FIG. 2: An example of a tt junction due to a rotating domain 
walls in a system with a d-wave order parameter. From left 
to right, the crystallographic axes rotates by 7r/6 across each 
domain wall. After three domain walls, the crystal axis has 
returned to itself, but the order parameter has changed sign. 



quently, the ground-state is a striped superconductor in 
which the period of the superconducting order parameter 
is twice the period of the lattice structure. 

This is a very contrived example, in which the origin of 
the PDW does not require any strong coupling effects - 
it is derived from the internal structure of the pair-wave 
function. For the point of principle, i.e. to establish the 
viability of a PDW as a phase of matter, this analysis 
is useful. However, for practical applications, a less con- 
trived mechanism is needed. 



IV. SOLVED MODELS 



elucidate the conditions under which this occurs. 

In the context of the cuprates, there have been several 
studies looking for a striped superconducting state in the 
t— J or Hubbard models. On the one hand, DMRG calcu- 
lations by White and Scalapino^ have consistently failed 
to find evidence in support of any sort of spontaneously 
occurring tt junctions. On the other hand, a number of 
variational Monte Carlo calculations have concluded that 
the striped superconductor is either the ground-state of 
such a model^, under appropriate circumstances, or at 
least close in energy to the true ground state— These 
variational calculations are certainly encouraging, in the 
sense that they suggest that there is no large energetic 
reason to rule out the existence of spontaneously occur- 
ring PDW order in strongly correlated electronic systems. 
However, since no such state has yet been observed in 
DMRG or other "unbiased" studies of these models, we 
believe the mechanism of formation of these states, and 
indeed whether they occur at all in physically reason- 
able microscopic models, remain unsolved problems. It 
is these problems that we aim to address. 



B. TT junctions from d-wave symmetry 

In the case in which the preferred superconducting 
state has a non-trivial internal pair-structure, as in a 
d-wave superconductor, it is well known both theoreti- 
cally and experimentally that tt junctions can be achieved 
by suitable orientation of the crystal fields across grain 
boundaries. This is not the physics we have in mind. 
However, since it is simple and well understood, we start 
with a highly artificial model in which d-wave symmetry 
gives rise to a striped superconductor. 

We consider the case of a d-wave superconductor in a 
square lattice, with a strong crystal field coupling which 
locks the lobes of the pair wave-function along the crys- 
tallographically defined x and y directions, as shown in 
Fig. O We imagine we have made a striped array in 
which the crystallographic axis rotates by 7r/6 between 
neighboring strips. Since rotation by 7r/2 is a symme- 
try, every third strip is identical. However, the d-wave 
order parameter changes sign under this rotation. Conse- 



A. Model of a single junction 

To start with, we consider a system consisting of three 
strips - a "left" superconductor described by a mean-field 
Hamiltonian with gap function AL(fc) (not to be con- 
fused with the order parameter of Eq. (|l.ip )and quasi- 



particle energy EL{k) ~ y [^i(^)]^ + I^L(fc)|^, a "right" 
superconductor which is also described by a mean-field 
Hamiltonian with the index L replaced by i?, and be- 
tween them a strongly correlated insulating " barrier" 
region. Here k is either a ID vector (if this is a line junc- 
tion) or a 2D vector if this is a planar junction. Note 
that we have assumed the left and right superconductors 
are thin in the transverse direction - otherwise we would 
have to include a transverse band index, as well. We as- 
sume that both superconductors, by themselves, preserve 
time-reversal symmetry, so in both cases the phase of the 
superconducting order (modulo tt) can be defined so that 
Aa(A:)e-*^° is real. 

The three decoupled subsystems are thus described by 
the Hamiltonian 



Hq — Hl + Hi 



(4.1) 



where, in the discussion below, we will consider several 
different examples for the insulating barrier. However, in 
all cases, we will assume that the superconducting Hamil- 
tonians have the quadratic forms 



ka 

+ V[A4fc)at at 



h.c. 



with a = L or R in which 



a r 



(4.2) 



(4.3) 



where cj^ - ^ creates an electron of spin polarization a at 
position f on, respectively, the left superconductor, the 
barrier, the right superconductor, for a = L, B, and R, 
and N = J2r^ ^he "area" of the junction. 
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The three subsystems are coupled together by a single- 
particle hopping term, 



H' 



r.a 



i'tLc[p,jCp^a + tRc'Rp^Cp^a + h.c] (4.4) 



(with the convention that cp^^ = CB,r,iy are the electron 
annihilation operators in the barrier region). Our pur- 
pose is to determine how the ground-state energy de- 
pends on the difference in phase between the supercon- 
ducting order parameter on the left and right side of the 
barrier, 9]^ — 0^. Specifically, we will consider the limit 
in which the matrix elements in H' tend to zero, so that 
the leading 9 dependent term in the ground-state energy 
can be computed in fourth order perturbation theory. 



A£; = -2N{tLtRf [Jcos{dL - Or) 



(4.5) 



where . . . as | | — ^ , and 2 1 t/fi l | ^ J is the Josephson 
coupling density across the barrier. 

An explicit expression can be obtained for J in terms 
of imaginary time-ordered correlation functions, as was 
shown in Ref . [3| ■ This can be derived either by making 
a Laplace transform of the ordinary perturbative expres- 
sion, or directly from an Euclidean path integral: 

J^^j dld2di d4f^i(l,2)J^^(4,3)f(l,2;3,4) 

(4.6) 

(4.7) 



where 1 = (ti, ri), 



dl 



E 



dri 



(in the limit [3 ^ oo) and 



i^a(l,2) 



Tr 



(4.8) 



r(l,2;4,3) = (T, ^J, _^(ri)4^^^(r2)c,-3,x(T3)c,^,,T(T4) 

It is convenient to express V as the sum of a non- 
interacting piece, corresponding to processes in which 
two electrons tunnel through the barrier independently of 
each other, plus a correction term, F, which expresses the 
effects of interactions within the harrier region between 
the tunneling electrons: 

f(l, 2; 3, 4) = Gt(1, 4)G|(2, 3) + r(l, 2; 3, 4) (4.9) 

where 

G.(l,2) ^ (r. [ct ,,(ri)c.^„.(r2)]) . (4.10) 
Finally, we define two contributions to J: 

J=Jl+J2 (4.11) 

where Ji is the non-interacting portion, and J2 is the 
portion proportional to F. 



B. General considerations 

Under most circumstances, Ji > since it is propor- 
tional to F^G^ and hence does not depend on the sign of 
G (see [45!). In Ref. '3, the situation in which there are 
attractive interactions between electrons in the barrier 
was investigated. In this case, J2 is also, generally, pos- 
itive. This earlier study dealt with the conditions under 
which J2 > Ji, leading to anomalously large values of the 
IcR product of the junction. In a sense, we are investigat- 
ing the converse problem, in which repulsive interactions 
lead to J2 < 0, while at the same time | J2I > Ji. 

In most cases with conventional superconductors, in 
which |A| is the smallest energy scale and ^0 (the super- 
conducting coherence length) is the largest length scale in 
the problem, it is typically the case that Ji | J2 1 . The 
dominant processes that contribute to Ji involve pairs 
tunneling within an (imaginary) time interval of each 
other ^ 1/|A| and within a radius ~ ^0 of each other. 
In contrast, the processes that dominate J2 involve elec- 
trons tunneling within an " interaction time" Tint and an 
interaction range, rint of each other. Thus, one generally 
expects I J2 / Ji I is small in proportion to a positive power 
of |A|ri„4 and rint/^o- In conventional superconductors, 
where these factors are very small, this is the deciding 
factor, so that except under very special circumstances, 
one expects positive ("ferromagnetic") coup ling across 
any junction. Indeed, as discussed in Ref. |4J, in the 
limit that A is small, Eq. (|4.6p reduces to the famil- 
iar form for a tunnelling Hamiltonian with an effective 
hopping <off = tbgiR across the barrier, where 



dl d2 

Nf3 



■G,(l,2). 



(4.12) 



Conversely, in "high temperature superconductors," in 
which A is not all that small, it is reasonable to expect, 
at least under some circumstances, that the coupling is 
negative - the long sought 7r-junction. 

There is one other factor that can suppress Ji relative 
to J2: When the two electrons tunnel across the barrier 
independently of each other, their individual crystal mo- 
menta are conserved in the process. Thus, if we remove a 
pair of electrons at momentum k and —k from the right 
superconductor, they must be injected into the left su- 
perconductor at the same momenta. If the right and left 
superconductors are the same, this is not a problem - 
electrons can be removed from near the Fermi energy on 
the left and added on the right at low energy. However, 
if the left and right superconductors have different values 
oi kp, then no such low energy process is possible. Con- 
versely, for correlated tunneling, only the center of mass 
momentum of the tunnelling electrons is conserved, so a 
pair of electrons removed from the left with momenta k 
and — fc, can be injected into the right superconductor 
with different momenta q and —q. 
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C. Explicit model "barriers" 

To make the considerations explicit, we will compute 
J for a model of the barrier in several limits. Specifically, 
we will consider the case in which the barrier consists of 
a Hubbard chain: 

Hb ^ ^ efi{r} - t ^ [4, o-Cf'.t^ + ti-' 

r <r,r'>,fT 

4,T4,iC-,iCr,T (4.13) 

r 

in the strong coupling limit, where the hopping matrix 
element t — > 0, with site energy e and on-site repulsion 
between two electrons on the same site, C/ > 0. We will 
consider two cases - the empty chain case in which e > 0, 
and the half-filled case, in which e<0<?7— |e|. In 
the half-filled case, the groundstate of the barrier is 2^ 
fold degenerate; we resolve this degeneracy by applying a 

staggered magnetic field, o- '^'''4o-'^Fdr where Q 

is the Neel ordering vector, and we take the limit h ~* 
at the end of the calculation. 

For simplicity, we will also assume that A l and A^^ are 
independent of /c, and that over the relevant range of en- 
ergies, the density of states in both superconductors can 
be approximated by a constant, pL and pjj, respectively. 



1. The empty barrier 

For the case of an empty barrier, e > 0, a straightfor- 
ward (but tedious) calculation reveals that 



changed. However, if the two superconductors are suffi- 
ciently different that iJ^, Efi ^ A, then 



Ji 




1 



1 



If the left and right superconductors are the same 



J I = 



P 
2eA 



f{y) - V 



(4.15) 



where 



/(ya) = 



cos ^(j/g) 



Va 



Ua = c/IAqI, and in this case y = vl ^ Vr- Here /(y) 
is a smooth function with /(O) = 7r/2, /(I) = 1, and 
fiy) ~ ln(2?/)/y as y -> oo. 

If the two superconductors have quite different values 
oi kp, then Ji gets substantially suppressed. Let El and 
Er be, respectively, the energy of the left quasiparticle 
at kp R and the energy of the right quasiparticle at fc^.i- 
So long as El, Er <C A, the magnitude of Ji is hardly 



Ji = PLfiVL) 



Eh 



PrUvr) 



(4.16) 



The interaction correction derives from terms that 
would have contributed to Ji, but are suppressed due 
to the repulsive interaction, and so has the opposite sign: 



J2 = -'2-f{yL)f{yR)pLPR 



U 



e(C/ + 2e) 



(4.17) 



In conventional superconductors, typically the band- 
width, W, is the largest energy in the problem, and hence 
p ^ is "small." Since Ji ^ p and J2 ~ p^, this 

means that Ji ^ | J2I. The exception to this rule occurs 
in cases in which AEp = mm[ER, El] is a substantial 
fraction of the bandwidth. 

However, in high temperature superconductors, we ex- 
pect to find that A is not too much smaller than W. In 
this case, it will depend on details whether or not Ji or 
J2 dominate. Clearly, this situation is more likely the 
larger U. In the case of tunnelling between two identical 
superconductors, and in the limit U 00, IJ2I > Ji so 
long as 



pA > Xc 



[f{y) - y] 
4/2(y)(i-2/^)' 



(4.18) 



For y ^ \, Xc ^ 2// In y, condition (|4.18p requires 
impossibly large gap scales. However, Xc = 19/96 for y ~ 
1 and Xc ^ l/27r as y ^ 0.~ Manifestly, Ji is reduced in 
magnitude and J2 is unaffected by a substantial value of 
AEp- For instance, if we consider the case in which the 
left and right superconductors differ only in the position 
of the Fermi surface, J2 dominates so long as 



(4.14) 



[pAEp]^f{y) > pA 



(4.19) 



2. The antiferromagnetic barrier 

We now turn to discuss the case e < < U ~ \e\, 
in which the barrier is half filled. As mentioned previ- 
ously, we introduce a staggered Zeeman field h to lift the 
groundstate degeneracy, and we take /i — > at the end 
of the calculation. The lowest order perturbative term in 
the Josephson coupling consists of all processes in which 
a pair of electrons tunnels through the barrier, subject to 
the constraint that the barrier has to return to its Neel 
ordered ground state at the end of the process. 

For simplicity, we consider the U —^ 00 case, in which 
no doubly occupied sites are allowed in the barrier. (We 
have obtained similar expressions for finite U, but they 
contain no qualitatively different physics.) In this case, 
processes in which electrons tunnel though different sites 
contribute to Ji and those in which they tunnel through 
the same site contribute to J2. As discussed in Ref. Is^, 
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serving term. (If particle- hole symmetry is present, this 
cancellation is exact.) 

We have evaluated the expressions in Eqs. (|4.20p and 
(j4.2ip numerically, for the case in which the two super- 
conducting strips are identical chains with a k indepen- 
dent gap function, Aa,k = A and the band-structure is 
the simplest tight-binding result, ea,k — — 2icos(fc) — ^. 
The results are summarized in Fig. [31 where the x axis is 
the average density of electrons (n) in the superconduc- 
tors, and the y-axis is the magnitude of A; the solid lines 
separate the region of net negative from the regions of net 
positive Josephson coupling. The general trends alluded 
to in the above discussion are clear from this figure. 



FIG. 3: Phase diagram for the antiferromagnetic barrier prob- 
lem, as a function of the pairing gap A (in units of the band- 
width W — 4:t) and the average density (n) in the supercon- 
ductors on either side of the barrier. When A is large or (n) 
is close to 1 (half filling), then a negative Josephson coupling 
(J < 0) is favored. 



the single site process necessarily involves an exchange 
of two fermions, and so makes a negative contribution to 
J. A new feature of the present problem is that there 
an antiferromagnetic "Umklapp" contribution to Ji, in 
which the tunneling electron exchanges momentum n/a 
with the antiferromagnet; this makes a negative contribu- 
tion to Ji which, under appropriate circumstances, can 
largely cancel the usual positive contribution. 

It is a straightforward exercise to obtain the resulting 
expressions: 



k,k' 
1 



^lS \ I ^R,k 



,'^^L,kJ \'^^B.,k 



_ ] ( ^k — k' ^ k—k' -\-7Ty 



EL,k + \^\ 



where Sj: is the Kronecker delta function, and 



2^^r k I \ '^Ej^^p 



(4.20) 



(4.21) 



Generically, the magnitude of the umklapp term (pro- 
portional to 5j:_j:,^^^) is considerably smaller than the 
momentum conserving term, since it is typically not pos- 
sible for both k or k + iry to be close to the Fermi surface 
unless a special "nesting" condition is satisfied. For ex- 
ample, if the two superconductors are one dimensional 
wires close to half filling, then the separation between 
the two Fermi points is indeed close to tt. In that case, 
the umklapp term tends to cancels the momentum con- 



D. A model with a striped superconducting 
ground state 

The basic ingredients necessary for finding a model 
with a striped superconducting ground-state are present 
in the above discussed examples. However, in the interest 
of making everything explicit, we consider the following 
problem. Firstly, we generalize the problem to c? + 1 di- 
mensions. We consider an array of parallel d dimensional 
hyper-planes with weak attractive interactions, and a 
particle-hole symmetric dispersion. Each hyperplane 
thus has a uniform superconducting order parameter. For 
weak attractions between electrons and large enough d, 
the superconducting correlations in each hyperplane can, 
presumably, be accurately treated in the context of BCS 
mean-field theory. The particle-hole symmetry implies 
that Er = Er , where n is the d-dimensional antiferro- 
magnetic Neel ordering vector. Sandwiched between each 
superconducting hyperplane there is a d-dimensional in- 
sulting hyperplane, which is the d-dimensional general- 
ization of the antiferromagnetic barrier discussed above. 
Finally, we couple neighboring hyperplanes with an arbi- 
trarily weak hopping matrix, t. The Josephson coupling 
per unit hyperarea is then computed as above. 

The particle-hole symmetry of the model insures that 
Ji = 0, and that therefore the Josephson coupling per 
unit hyper-area is proportional to J2, and hence is neg- 
ative. Thus, the ground-state of this system is a state 
with coexisting SDW and striped superconducting or- 
der. Note that the phase of the antiferromagnetic order 
is determined by the sign of h. If we take ft. > in 
all cases, then the antiferromagnetic order is "in-phase." 
In this case, the Hamiltonian is invariant under transla- 
tions in the direction perpendicular to the hyperplanes by 
the distance between insulating hyperplanes, and hence 
the striped superconductor spontaneously breaks trans- 
lational symmetry. However, we can also imagine stag- 
gering the sign of h from one hyperplane to the next, so 
as to mimic the antiphase spin-stripe order seen in the 
cuprates. In this case, the underlying unit cell is twice as 
large, and the striped superconductor does not sponta- 
neously break translational symmetry. However, in both 
cases, the spatial average of any component of the super- 
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conducting order parameter is 0, so this state conforms 
to our definition of a PDW. 

This model is, admittedly, somewhat contrived, espe- 
cially in the assumed particle-hole symmetry. However, 
since at this point J2 has a finite magnitude, the state is 
robust with respect to small deformations of the model. 

The construction of solvable microscopic models with 
a striped superconducting ground-state, and the elucida- 
tion of the ingredients necessary to get anti-phase super- 
conducting order across a correlated insulating barrier, 
are the principle results of the present paper. 



E. Quasiparticle spectrum of a striped 
superconductor 

The spectrum of an s-wave superconductor is typically 
fully gapped, while a d-wave superconductor has a gap 
that vanishes only on a few isolated points (lines in 3d) 
on the Fermi surface. These features persist even in the 
presence of coexisting CDW or SDW order of various 
sortsl-. In particular, the density of quasiparticle states 
is either zero for energies less than a non-zero minimum 
value (in the s-wave case), or vanishes (linearly) as the 
energy approaches the Fermi surface. This effect of su- 
perconducting order is in contrast to most other orders, 
including SDW and CDW order, in which, so long as 
the effective ordering potential is weak compared to the 
Fermi energy, at most a portion of the Fermi surface is 
gapped leaving the system with a (smaller) reconstructed 
Fermi surface, but a Fermi surface none-the-less. 

The gap character in a PDW is more similar to that 
of a CDW than to a uniform superconductor. Specif- 
ically, so long as the order parameter is not too large, 
an ungapped reconstructed Fermi surface remain o^^'^^ . 
Thus, a PDW typically has a finite density of states in 
the superconducting phase. 

For band-structure parameters characteristic of a typ- 
ical underdoped or optimally doped curpate, the quasi- 
particle spectrum of a striped superconductor has a large 
(d- wave-like) gap near the " antinodal points" {i.e. near 
(tt, 0) and (0, tt)), but has a closed, Fermi-surface pocket 
which closely follows the contours of the portion of the 
normal-state Fermi surface near the " nodal points" [i.e. 
where the Fermi surface crosses the (0, 0) to (tt, tt) cord]. 
In short, it closely resembles certain phenomenological 
descriptions'*'' of the "pseudo-gap." The low energy spec- 
tral weight is substantial only on a portion of the bare 
Fermi surface near the nodal (tt, tt) direction^i, forming 
a "Fermi arc"i^. 



V. COUPLED ORDER PARAMETERS 

In this section, we explore the aspects of the theory of 
a PDW that can be analyzed without reference to mi- 
croscopic mechanisms. We focus on the properties of or- 
dered states at T = 0, far from the point of any quantum 



phase transition, where for the most part fluctuation ef- 
fects can be neglected. (The one exception to the general 
rule is that, where we discuss effects of disorder, we will 
encounter various spin-glass related phases where fluctu- 
ation effects, even at T = 0, can qualitatively alter the 
phases.) For simplicity, most of our discussion is couched 
in terms of a Landau theory, in which the effective free 
energy is expanded in powers of the order parameters; 
this is formally not justified deep in an ordered phase, 
but it is a convenient way to exhibit the consequences of 
the order parameter symmetries. 

Specifically, the emphasis in this section is on the inter- 
relation between striped superconducting order and other 
orders. There is a necessary relation between this order 
and CDW and nematic (or orthorhombic) order, since the 
striped superconductor breaks both translational and ro- 
tational symmetries of the crystal. From the microscopic 
considerations, above, and from the phenomenology of 
the cuprates, we also are interested in the relation of su- 
perconducting and SDW order. We assume that the host 
crystal is tetragonal, and that there are therefore two po- 
tential symmetry related ordering wave-vectors, Q and 
Q, which are mutually orthogonal. Spin-orbit coupling 
is assumed to be negligible. We introduce the following 
order parameter fields: the nematic, N , the PDW Ag, 
the SDW Sq, and the CDW pk (where K ^ 2Q), and, 
of course, the corresponding orders at the symmetry re- 
lated wave-vector, Q. The Landau effective free energy 
density can then be expanded in powers of these fields: 

T = T2+:F3 + T4 + ... (5.1) 

where the quadratic term, is simply a sum of decou- 
pled terms for each order parameter, 

•^3 = 1s[p*kSq ■ Sq + pISq • Sq + C.C.] (5.2) 
+7aKA* qAq + p:^A* qAq + c.c] 
+5AiV[A^AQ + A* qA_q - A^Aq - A* ^A_q] 

+9cN[p*jiPK - P*kPk]^ 

and the fourth order term, which is more or less standard, 
is shown explicitly in Appendix [XI (Detailed discussions 
of the microscopic definition of the PDW order parameter 
is also discussed in this Appendix.) 

The effect of the cubic term proportional to 7^ on 
the interplay between the spin and charge components 
of stripe order has been analyzed in depth in Ref. [sol . 
Similar analysis can be applied to the other terms. In 
particular, it follows that the existence of superconduct- 
ing stripe order (Ag ^ 0, and Ag — 0), implies the exis- 
tence of nematic order {N ^ 0) and charge stripe order 
with half the period {p2Q 7^ 0). However, the converse 
statement is not true: while CDW order with ordering 
wave-vector 2Q or nematic order tend to promote PDW 
order, depending on the magnitude of the quadratic term 
in JF2, PDW order may or may not occur. 
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One new feature of the coupling between the PDW 
and CDW order is that it produces a sensitivity to disor- 
der which is not normally a feature of the superconduct- 
ing state. In the presence of quenched disorder, there 
is always some amount of spatial variation of the charge 
density, p{r) , of which the important portion for our pur- 
poses can be thought of as being a pinned CDW, that is, 
a CDW with a phase which is a pinned, slowly varying 
function of position, p{r) = \pk \ cos[K ■ r + (j){r)]. Be- 
low the nominal striped superconducting ordering tem- 
perature, we can similarly express the PDW order in 
terms of a slowly varying superconducting phase, A(r) — 
\AQ\exp[iQ -r + ieQir)] + \A^Q\exp[-iQ ■ r + ie^gir)]. 
The resulting contribution to J-3 is 

^3,7 = 27a|pj^AqA_q| cos[2e_{r) - </)(r)]. (5.3) 

where 

0±ir) ^ [0Qir)±e^Qir)]/2; (5.4) 
0±Q{r) = [0+(r)±0_(r)]. 

The aspect of this equation that is notable is that the dis- 
order couples directly to a piece of the superconducting 
phase, 9- . No such coupling occurs in usual momentum 
superconductors. 

It is important to note that the condition that A(r) be 
single valued imposes a non-trivial topological constraint 
on possible vortices. Specifically, an isolated half-vortex 
(about which the phase winds by tt) is forbidden in cither 
6'-|_ or 6'_; vortices must occur either as full (27r)-vorticcs 
in one or the other phase field, or as a bound pair of a 0+ 
and a 0_ half-vortex. An important consequence of this 
phase coupling is that the effect of quenched disorder, 
as in the case of the CDW itself, destroys long-range 
superconducting stripe order. (This statement is true^, 
even for weak disorder, in dimensions d < 4.) Naturally, 
the way in which this plays out depends on the way in 
which the CDW state is disordered. 

In the most straightforward case, the CDW order is 
punctuated by random, pinned dislocations, i.e. 2t: vor- 
tices of the (f) field. The existence of the coupling in 
Eq. (|5.3p implies that there must be an accompanying 
TT vortex in The condition of single- valued-ness im- 
plies that there must also be an associated half-vortex 
or anti-vortex in 9+—. If these latter vortices are fluc- 
tuating, they destroy the superconducting state entirely, 
leading to a resistive state with short-ranged striped su- 
perconducting correlations. If they are frozen, the re- 
sulting state is analogous to the ordered phase of an XY 
spin-glass: such a state has a non-vanishing Edwards- 
Anderson order parameter, spontaneously breaks time- 
reversal symmetry, and, presumably, has vanishing re- 
sistance but no Meissner effect and a vanishing critical 
current. In 2D, according to conventional wisdom, a spin- 
glass phase can only occur at T = 0, but in 3D there can 
be a finite temperature glass transition^'^. 

In 3D there is also the exotic possibility that, for weak 
enough quenched disorder, the CDW forms a Bragg-glass 



phase, in which long-range order is destroyed, but no free 
dislocations occur^^ i^^i^^ . In this case, (j) can be treated 
as a random, but single-valued function - correspond- 
ingly, so is 9-. The result is a superconducting Bragg- 
glass phase which preserves time reversal symmetry and, 
presumably, acts more or less the same as a usual super- 
conducting phase. It is believed that a Bragg-glass phase 
is not possible in 2D^. 

A summary of the characterization of the PDW, SC 
and CDW phases in terms of their order parameters, their 
inter-relations and their sensitivity to quenched disorder 
appears in Table |T1 

Another perspective on the nature of the supercon- 
ducting state can be obtained by considering a compos- 
ite order parameter which is proportional to AqA_q. 
There is a cubic term which couples a uniform, charge 
4e superconducting order parameter, A4, to the PDW 
order: 

J-^ = .g4{A^[AQA_Q + AqA_q] + c.c.} (5.5) 

This term implies that whenever there is PDW order, 
there is also necessarily charge 4e uniform superconduct- 
ing order. However, since this term is independent of 
it would be totally unaffected by Bragg-glass formation 
of the CDW. The half-vortices in 9+ discussed above can 
simply be viewed as the fundamental (hc/Ae) vortices of 
a charge 4e superconductor. 

Turning now to the quartic terms in ^4, there are a 
number of features of the ordered phases which depend 
qualitatively on the sign of various couplings. Again, 
this is very similar to what happens in the case of CDW 
order - see, for example, Ref. [13] ■ For instance, de- 
pending on the sign of a certain biquadratic term, either 
unidirectional (superconducting stripe) or bidirectional 
(superconducting checkerboard) order is favored. 

Finally, we comment on the case of coexisting uniform 
and striped superconducting order parameters. Such a 
state is not thermodynamically distinct from a regular 
(uniform) superconductor coexisting with a charge den- 
sity wave, even if the uniform superconducting compo- 
nent is in fact weaker than the striped component. There- 
fore, we expect many of the special features of the striped 
superconductor (such as its sensitivity to potential dis- 
order) to be lost. In Appendix [BJ we extend the Landau 
free energy to include a uniform superconducting com- 
ponent, and show that this is indeed the case. 

VI. MODEL FOR A LAYERED STRIPED 
SUPERCONDUCTOR 

We now discuss a low energy effective theory for a lay- 
ered, disordered striped superconductor. As we saw in 
the previous section, disorder inevitably nucleates half 
vortices which are pinned to dislocations in the CDW. 
The system therefore has Ising-like degrees of freedom 
which are the "charges" (or vorticities) of these tt vor- 
tices. We assume that their positions are quenched, ran- 
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dom variables. As we discussed in Ref. [2^, a PDW in 
a layered superconductor can give rise to a frustration 
of the inter-layer Josephson coupling, given a structure 
in which the CDW is rotated by 90° between layers [as 
occurs in the low temperature tetragonal (LTT) phase of 
La2-2;Ba2;Cu04]. To make the model relevant to this sys- 
tem, we consider the case where the inter-layer Josephson 
coupling vanishes identically. 

Under these assumptions, the interaction energy be- 
tween half vortices is composed of two parts: the mag- 
netic energy and the kinetic energy associated with the 
screening currents that surround the vortices. The to- 
tal interaction energy for a given configuration of half 
vortices is^ 



ijaf3 a(5 



(6.1) 

Here, qia — ±1 and Ria are the "charge" (or vorticity) 
and the position of the ith vortex in layer a, respectively, 
Za is the z coordinate of layer a, and Qa = J2i lia- The 
interaction potentials u{R) and v (z) are given by 



u R] = 



Mo 



(2^ 



2A a (k) fc2 



Jk-p 



1 



where p and z are the radial and the z axis separation be- 
tween the two vortices, respectively, a (k) 
Bo 



1 



(Afc)^ 



4eA 



, Saf3 is a Kronecker delta function of the layer 



indices a, (3 of the two vortices, A is the in-plane pene- 
tration depth, d is the inter-layer distance, and 
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where L is a long wavelength cutoff of the order of the 
linear dimension of the system, ^ is the coherence length. 
El {x) = JJj" ^~^dt is the exponential integral function, 
and 7 is Euler's constant. The second term in Eq. (|6.ip 
diverges as InL, unless Qa — Q — const, for all a. 
In the absence of an external magnetic field, this term 
constrains Qa — 0. This is the usual "charge neutrality" 
condition which comes from the infrared divergence of 
the vortex self energy. 

For \J ^ z"^ :$> A, Eq. (|6.2p reduces to 



(r).^ 



Sap - ^e- 



mi. 

p 



(6.4) 



The statistical mechanics problem of a finite density of 
half vortices with quenched random positions, and whose 
interaction is given by Eq. (|6.ip . is an interesting un- 
solved statistical mechanics problem.— 

As we mentioned, the putative superconducting glass 
phase (in which the half vortices are frozen) necessar- 
ily involves time reversal symmetry breaking. It should 



therefore be detectable by measuring the magnetic fields 
associated with the spontaneous half vortices which oc- 
cur at dislocations in the charge order. For example, 
consider a half vortex at the surface of the sample (taken 
to be at z = 0). The asymptotic form of the magnetic 
field above the surface in the limit < z <C A, is 



(z « A) 



■d 

8A2e I z 



(6.5) 



where d is the interlayer distance. [For z > A, screening 
by diamagnetic currents in the other planes becomes con- 
siderable, and crosses over to ~ (p')-] Assuming 

A « 2000 A and d « 15 A, Eq. (gSl) gives B^ w 300 

G. At a distance of about 1000 A from surface, the result- 
ing characteristic fields of the order of 0.3 G (assuming a 
half vortex right at the surface) are well within the reso- 
lution of current local magnetic field measurement tech- 
niques. The onset of a spontaneous, random magnetic 
fields of this order at the glass transition temperature 
(which is presumably also signaled by the vanishing of 
the linear resistivity) would be a dramatic confirmation 
of the PDW scenario in \ doped LBCO. 



VII. FINAL THOUGHTS 

The strong suppression of the three-dimensional super- 
conducting Tc in La2-xBaa;Cu04 at a; = 1/8, and other 
"1/8 anomalies" in the La2-xSrxCu04 family of cuprate 
superconductors, have long been interpreted in the liter- 
ature as evidence that charge order competes with high 
temperature superconductivity. (More recentl y^^'^^'^'^ , 
clear evidence of a 1/8 anomaly has been adduced in 
YBa2Cu307_5, as well.) However, the remarkable prop- 
erties of La2_a;Baa;Cu04, particularly the fact that the 
anti- nodal gap is largest at x — 1/8, where the dynam- 
ical layer decoupling is observed, strongly suggests that 
charge stripe order can he part of the mechanism of su- 
perconductivity as argued in Ref. [6^ . 

The superconducting state presented here represents 
a new face of the interplay between superconductivity 
and spin and charge order. If the SC and SDW orders 
simply competed, the PDW would be a rather unnatural 
state; it is natural if self-organized inhomogeneities are 
an essential feature of the mechanism. In this sense, the 
PDW state should also be regarded as an electronic liquid 
crystal phase. ^'^ Numerous and important implications 
follow from the properties of this state. 



A. Implications 

Glassy superconductors: The most dramatic conse- 
quence of the PDW physics is that, in the presence 
of weak, quenched disorder, the superconducting phase 
gives way to a regime of glassy behavior, where strong 
local superconducting correlations extend up to a finite 
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correlation length set by the strength of the disorder. In 
the regime of temperature below the onset of substantial 
local superconducting coherence but above the transition 
to a fully superconducting state (if one occurs), the sys- 
tem can be characterized as a "failed superconductor" 
Clearly, in this regime, the longitudinal resistivity will be 
small compared to normal state values, as will the quasi- 
particle contributions to the thermopower and linear Hall 
resistance. Moreover, one expects to see strong local indi- 
cations of superconductivity, especially the formation of 
a superconducting pseudo-gap in the single-particle spec- 
trum with a form similar to that of the ordered phase. 

Taking into account that the materials are ultimately 
3D, there will generally be a glass transition at T = 
in this regime. In the absence of an applied magnetic 
field, the glass phase can be characterized by its spon- 
taneously broken time-reversal symmetry. Presumably, 
the glass phase has zero resistance, zero critical current, 
and fails to exhibit a full Meissner effect. More gener- 
ally, for a range of temperatures, including temperatures 
above Tg, the magnetic response of the PDW will be 
characterized by slow dynamics, a broad distribution of 
relaxation rates, and probably a certain degree of history 
dependence; all the dramatic and confusing features of 
spin-glasses, but amplified by the coherent orbital orga- 
nization of a superconducting state. 

Interlayer decoupling: When a PDW state occurs in a 
quasi-2D (layered) material, it will frequently be the case 
(depending on the interlayer geometry) that the usual 
Josephson coupling between neighboring layers vanishes. 
This feature was explored at length in our earlier pa- 
per, Ref. [l^]. In the case of the LTT structure of 
La2-a;Ba2;Cu04, the fact that the stripes in neighbor- 
ing planes run at right angles to each other insures that 
(in the absence of disorder), the Josephson coupling be- 
tween neighboring planes vanishes. Strictly speaking, 
this does not mean that there is no coupling between 
plains, as would occur in a putative "fioating phase" 
There are always couplings between farther neighbor lay- 
ers, and higher order Josephson couplings between neigh- 
boring planes (which couple A4 in 3D), and these inter- 
actions are relevant in the renormalization group sense at 
any temperature below the putative Kosterlitz-Thouless 
(KT) transition temperature. However, at the very least, 
it means that the 3D superconducting state is enormously 
more anisotropic than would be expected on the basis of 
the bare electronic anisotropy. Moreover, as long as su- 
perconducting coherence in a given plane is limited due 
to quenched disorder, if these higher order couplings are 
weak enough, they can be essentially ignored. 

Striped superconductors in strongly correlated models: 
While we have established, as a point of principle, that 
well defined models exist that support a PDW phase, it 
still makes sense to make predictions that can be tested 
in "numerical experiments" on i — J and Hubbard mod- 
els. The fact that spontaneously occurring 7r-junctions 
have not yet been reported in extensive previous DMRG 
studies^ is worrisome, in this regard. 



Based upon the insights gained in the above, we pro- 
pose DMRG studies of microscopic models in which we 
expect indications of PDW formation can be observed. 
Two examples are: 

i) A three leg Hubbard ladder in which the outer two 
legs have a negative U and the inner-leg a positive U . 
The chemical potential and on-site energies should be 
chosen so that the inner leg is near half-filling, and so 
strongly antiferromagnetically correlated. The density of 
electrons on the outer ("superconducting") legs can be 
varied, and need not be identical. This model is thus a 
close relative of the model that we treated in Sec. VI. C. 2, 
except that we do not treat the superconducting legs in 
mean-field theory and we are not restricted to consider- 
ing parametrically small coupling between the legs. The 
tendency to a striped SC phase can be tested by studying 
the sign of the pair-field correlations between the upper 
and lower legs. When the outer legs are near half filling, 
we expect to see negative correlations, indicative of PDW 
formation. Indeed, we have already seen such behavior 
in preliminary DMRG studies of this modelfi-^'. 

ii) A five leg Hubbard ladder with all repulsive interac- 
tions, constructed to consist of two outer two-leg ladders 
and an inner Hubbard chain. Again, the chemical po- 
tential and on-site energies should be chosen so that the 
inner chain is near half-filling, making it strongly antifer- 
romagnetically correlated, while the density of electrons 
on the outer ladders can be varied, and need not be iden- 
tical. By making the densities on the outer legs different 
enough, we are confident that a PDW state can be in- 
duced. However, we are very curious to learn how robust 
this state is under less contrived conditions. We are now 
undertaking such an investigation^. 



B. Speculations 

Striped SC phases in La2-xBaxCuOi at x = 1/8: We 
were originally motivated in this study by experiments 
in La2-a:Ba2;Cu04, and we remain optimistic that they 
are, indeed, evidence of the existence of a PDW phase. 
In order to make more than an impressionistic compar- 
ison with experiment, the nature of the superconduct- 
ing glass phase will need to be understood theoretically, 
much better than we do at present. However, in the ab- 
sence of such theoretical control, we can still make some 
speculative statements concerning the relation between 
experiment and theory. 

Tentatively, we would like to identify the point at 
which resistivity vanishes in the c-direction as the true 
point of the 3D glass transition, Tg « lOK. If this is the 
case, the linear-response resistivity at lower temperatures 
is truly zero, although various non-linear and hysteretic 
processes may complicate the measurements. We expect 
very long time-scales and a degree of history dependence 
of macroscopic measurements to begin being significant 
at considerably higher temperatures, T > Tg. 

It is natural to associate the point at which the in- 
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plane resistivity apparently vanishes Tkt ~ 16K, with 
what would have been the KT transition of a single plane 
in the absence of quenched disorder. We do not think 
this is a true transition, and would expect that if experi- 
ments could be carried out with higher precision than has 
currently been possible, the in-plane resistivity would be 
found to have a finite value for Tg < T < Tkt- In 
relatively clean samples, simple scaling arguments sug- 
gest that the residual resistivity should be proportional 
to where ^ is the coherence length, which in turn is 
roughly the distance between dislocations. 

The sharp drop in the resistivity at the spin-ordering 
temperature, Tsdw ~ 42K, is probably not a true phase 
transition, either, but rather marks the sudden onset of 
significant intermediate scale superconducting coherence. 
However, we suspect that local pairing correlations per- 
sist to higher temperatures; as long as the spin order is 
strongly fluctuating at T > Tsdw, we imagine that phase 
coherence between neighboring superconducting stripes 
is prevented, while for T < Tsdw, antiphase supercon- 
ducting correlations extend over multiple stripes. We fa- 
vor this viewpoint for several reasons, most importantly 
because the gap-features seen in ARPES persist to higher 
temperatures f22^ 

ARPES spectrum of La2-xBaxCu04: A recent studyS'^ 
of the ARPES spectrum of La2-xBaa;Cu04 confirmed the 
conclusion of an earlier^^ ARPES/STM study concerning 
the existence of a gap consistent with a generally d-wave 
angle-dependence. However, the higher resolution study 
revealed that this gap has what appears to be a two- 
component structure, in that the gap near the anti-nodal 
point at (tTjO), is considerably larger (by a factor of 2- 
3) than a simple extrapolation from the nodal direction 
would suggest. Since there are several forms of density 
wave order known to be present in La2-a;Baa;Cu04, it is 
not straightforward to unambiguously identify particular 
spectral features with particular types of order. This 
is particularly problematic, since fluctuating order can 
also lead to a pseudo-gap with many similarities to the 
gap that would be produced in the corresponding ordered 
phase. 

That having been said, it is striking how much the ob- 
served spectrum in La2-xBaa;Cu04 resembles the mean- 
field spectrum found under the assumption that there 
is large amplitude PDW (striped) order coexisting with 
small amplitude uniform superconducting order. The 
PDW has a large gap along the anti-nodal [(7r,0) or (0,7r)] 
direction, and a Fermi pocket in the nodal [(7r,7r)] direc- 
tion, whose spectral weight is considerable only along an 
open "Fermi arc" region which nearly coincides with the 
bare Fermi surface. If, in addition, there is a uniform d- 
wave component to the order parameter, this completes 
the gapping of the Fermi surface. A spectrum that resem- 
bles the measured ARPES spectrum in La2-a;Baj;Cu04 
can be obtained^^ under the assumption that the PDW 
has a gap magnitude ApDw = 40meV, and uniform gap 
magnitude = 8meV. The ARPES spectrum is always 
measured at temperatures well above the bulk supercon- 



ducting Tc = so to the extent that this identification 
is correct, what is being observed is a pseudo-gap. More- 
over, as mentioned above, the gap structure persists to 
temperatures above Tsdw, although at the higher tem- 
peratures, the gap in the nodal regions decreases ap- 
preciably. Much more detailed analysis of the energy, 
momentum, and temperature dependence of the ARPES 
spectrum will be necessary to test this interpretation. 

Relevance to other cuprates: We have summarized 
the evidence for a state with the symmetries of a 
PDW state (or the glassy version of it) in i doped 
La2-2;Baa:Cu04. There is also evidence for dynamical 
layer decouphng effects in Lai.6-a;Ndo.4Sra;Cu04^, and 
also in La2-sSr2,Cu04 in a magnetic field i^^'^^ (For the 
possible relevance of these ideas to certain heavy fermion 
and organic superconductors see [TQI-) 

The quasiparticle spectrum of a PDW state has strik- 
ing qualitative resemblance to the spectra seen in ARPES 
in other cuprates, especially BSCCO 2212 and 2201. This 
possibly sheds new light on the issue of the "nodal-anti 
nodal dichotomy" : According to this interpretation, both 
the nodal and the anti-nodal gaps are superconducting 
gaps, the first being uniform and the other modulated. 
There are two distinct types of order ("two gaps"), but 
they are both superconducting, and so they can smoothly 
evolve into one another. Perhaps, as the temperature is 
decreased, the PDW gradually decreases and the uniform 
order parameter increases, while their sum (which is de- 
termined by relatively high energy microscopic physics) 
remains approximately constant. (Note that an early 
study^** of modulated structures seen in STM^^ii^ con- 
cluded that they could be understood in terms of just 
such a two-superconducting-gap state.) 

More generally, one of the most remarkable features 
of the pseudo-gap phenomena is the existence of effects 
of superconducting fluctuations, detectable2Ii2^i2S for in- 
stance in the Nernst signal, over a surprisingly broad 
range of temperatures and doping concentrations. At a 
broad-brush level^, these phenomena are a consequence 
of a phase stiffness scale that is small compared to the 
pairing scale. However, it is generally difficult to under- 
stand the existence of such a broad fluctuational regime 
on the basis of any sensible microscopic considerations. 
The glassy nature of the ordering phenomena in a PDW 
may hold the key to this central paradox of HTC phe- 
nomenology, as it gives rise to an intrinsically broad 
regime in which superconducting correlations extend over 
large, but not inflnite distances. 
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APPENDIX A: ORDER PARAMETERS 

We will now define the various order parameters intro- 
duced in Sec. |V] and discuss their symmetry properties. 
The striped superconducting order parameter Aq is a 
charge 2e complex scalar field, carrying momentum Q. 
To define it microscopically, we write the superconduct- 
ing order parameter as 



A(r,r') ^ (V^|(r)^|(f')) 



Ao + Age' 



iQ-R 



(Al) 



where R = {f+r')/2, F (f — r') is some short range 
function (for a "d-wave like" striped superconductor, 
F (r) changes sign under 90° rotation), and Aq is a uni- 
form order parameter. In the remaining of this appendix, 
we set Aq — 0. The effect of Aq is discussed in Appendix 
IbI In cases where a Aq is also non-zero (where Q is re- 
lated to Q by a 90° rotation) analogous terms with Aq 
replaced by Aq have to be added. 

The order parameters that may couple to Aq and their 
symmetry properties are as follows. The nematic order 
parameter TV is a real pseudo-scalar field; the CDW pK 
with K = 2Q is a complex scalar field; Sq is a neu- 
tral spin- vector complex field. All these order parameters 
are electrically neutral. Under spatial rotation by 7r/2, 
N -N, PK PR, Sq %, and Aq -> ±Aq, where 
± refers to a d-wave or s-wave version of the striped su- 
perconductor. Under spatial translation by r, iV — > iV, 
PK ^ e'^ Vif, Sq ^ e'Q-^SK, and Aq ^ c^^'^Aq. 

With these considerations, we write all the possible 
fourth order terms consistent with symmetry: 



•^4 = +u (5q ■ SqA*qA^q+Sq ■ SqA*qA_q + c.c.) 

+ [v+[S*Q -Sq + S^q- Sq] + i+[\pK\' + \pk\']) (|AqP + |A_Qp + |AqP + |A_Qp} 
-f (y^S^Q -Sq-S^q- Sq] + v^IpkI' - Ipa-P ]) (|AqP + |A_qP - |AqP - |A_qP) 
+.A^^ {(|AqP + IA^qH + (|AqP |A_Qp)} 

+A+ {(IAqI^ + \A.QfY + (|Aq|2 + |A_q|^)^} + A_ {(|Aq|^ - |A.qH^ + (|Aq|^ - |A-q|')'} + ■ ■ ■ 

(A2) 



where we have explicitly shown all the terms involving 
Aq, while the terms . . . represent the remaining quartic 
terms all of which, with the exception of those involving 
A'', are exhibited explicitly in Ref. Isol 

On physical grounds, we have some information con- 
cerning the sign of various terms in J- 4. The term propor- 
tional to u determines the relative phase of the spin and 
superconducting stripe order - we believe u > which 
thus favors a 7r/2 phase shift between the SDW and the 
striped superconducting order, i.e. the peak of the su- 
perconducting order occurs where the spin order passes 
through zero. The other interesting thing about this term 
is that it implies an effective cooperativity between spin 
and striped superconducting order. The net effect, i.e. 
whether spin and striped superconducting order cooper- 
ate or fight, is determined by the sign of |m| — w+ — u-, 
such that they cooperate ii \u\ — v+ — V- > and oppose 
each other if |u| — 11+ — < 0. It is an interesting pos- 
sibility that spin order and superconducting stripe order 



can actually favor each other even with all " repulsive" 
interactions. The term proportional to A_ determines 
whether the superconducting stripe order tends to be 
real (A_ > 0), with a superconducting order that sim- 
ply changes sign as a function of position, or a complex 
spiral, which supports ground-state currents (A_ < 0). 



APPENDIX B: COEXISTING UNIFORM AND 
STRIPED ORDER PARAMETERS 

In this appendix, we analyze the coupling of a striped 
superconducting order parameter Aq to a uniform order 
parameter, Aq. In this case, we have to consider in ad- 
dition to the order parameters introduced in Sec. |V] a 
CDW order parameter with wavevector Q, denoted by 
Pq. The additional cubic terms in the Ginzburg-Landau 
free energy are 
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.9p 



■k + c.c. 



(Bl) 



Eq. (|Bip shows that if both Aq and Ag are non-zero, there is necessarily a coexisting non-zero pg, through the 7q 
term. The additional quartic terms involving Aq are 



^4,n = MA(A^2AQA_Q + A*2AgA_Q + c.c.) +5|AoH|Aq|2 + |Aq| 



I An 



pq\ 



P2q| 



\p-2q\ 



(B2) 



Let us now consider the effect of quenched disor- 
der. Following the discussion preceding Eq. (|5.4|) . we 
write the order parameters in real space as A (r) = 
|Ao| e*®»-h|AQ| e^t^Q+Q-^' + IA.Ql e^(«-Q-e-'-) and p(r) = 
\pq\ cos {Q - r + (j)Q) -I- \p2q\ cos {2Q ■ r + (f>). Let us as- 
sume that the disorder nucleates a point defect in the 
CDW, which in this case corresponds to a 27r vortex in 
the phase (/)q. By the gp term in Eq. (jBip . this in- 
duces a 47r vortex in 0. (Note that a in the presence 



oi Pq, a 27T vortex in is not possible.) The 7a term 
in Eq. ()5.3p then dictates a 2tt vortex in the phase 
^- = (^Q ^ ^-q) /2- However, unlike before, this vor- 
tex does not couple to the global superconducting phase 
6+ ~ {9q + 9-q) /2. Therefore, an arbitrarily small uni- 
form superconducting component is sufficient to remove 
the sensitivity of a striped superconductor to disorder, 
and the system is expected to behave more or less like a 
regular (uniform) superconductor. 
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